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1 Introduction 

Understanding the structure of hadrons and their 
excitations stiU remains as an outstanding and fun- 
damental challenge in the hadron physics. The ex- 
istence of the hadrons is a direct consequence of the 
confinement of quarks and gluons within the quan- 
tum chromodynamics (QCD). Thus exploring the 
hadron structure is nothing but making clear the non- 
perturbative nature of QCD. 

Recent experimental activities at electron-beam 
facilities such as JLab, Bonn, Mainz, SPring-8, and 
GRAAL have brought new light to the study of ex- 
cited nucleon (N*) states^. Those facilities provide 
a considerable amount of precise data of the meson 
photo- and electro-production reactions on the nu- 
cleon target, and open a great opportunity to make 
quantitative investigations of the N* structure. 

Against the background of those experimental 
progress, continuous effort to extract the N* prop- 
erties from the world data of meson production ttN, 
"fN, and N{e,e') reactions is being made in Excited 
Baryon Analysis Center (EBAC) at Jefferson Lab. 
The analysis is pursued with a dynamical coupled- 
channels (DCC) approach proposed in Ref. The 
approach treats most relevant unitary cuts below 
W = 2 GeV, i.e., ttN, rjN, and ttttN which has un- 



stable ttA, pN, aN channels, in solving the coupled- 
channels equations. 

The objective of EBAC is more than just perform- 
ing the partial- wave analysis of the meson production 
reactions. We not only try to extract the N* param- 
eters, but also to map out the quark-giuon substruc- 
ture of the N* states. It will require comprehensive 
study combined with various hadron structure cal- 
culations such as constituent quark models, covari- 
ant models based on Dyson-Schwinger equations, and 
Lattice QCD simulations. The N* information ex- 
tracted from the reaction model analysis, such as that 
performed in EBAC, is vital for bridging the gap be- 
tween the actual reaction data and the hadron struc- 
ture calculations. 

The main subjects in EBAC are summarized as 
follows: 

1. Establish baryon spectrum. 

2. Extract N* parameters, in particular the elec- 
tromagnetic N-N* transition form factors, from 
analyzing the world data of meson production 
ttN, jN, N{e,e') reactions. 

3. Develop a method to connect the extracted 
form factors to the hadron structure calcula- 
tions and deduce the structure of the N* states. 
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In this paper we present an overview of current 
status and future plans of the EBAC program. In 
Sec. [2] we briefly explain a framework of the EBAC- 
DCC model, and present the results of our analysis of 
the meson production reactions in Sec. [3l In Sec. [4] 
we present the extracted N* pole positions within the 
current EBAC-DCC model. Finally we discuss our 
future plan in Sec. [5j 

2 EBAC-DCC model 

The EBAC-DCC analysis is based on a multi- 
channels and multi-resonances model'^' within which 
the partial wave amplitudes of + B{—p) 

M'{p') + B'{—p') are calculated from the following 
coupled-channels integral equation (suppressing the 
angular momentum and isospin indices): 

Tmb,ai'b'{p,p';E) = Vmb,m'b'{p,p';E) + 

dqq^VMBM" B"{p,q;E) X 

Gm-b" (<?; E)Tm"b",m'b' {q,p';E). (1) 

Here Vmb.m'b' and Gmb are the MB M'B' transi- 
tion potential and the MB Green function described 
below; the M"B" summation represents the coupled- 
channels effects among the ttN, rjN, ttA,pN, aN chan- 
nels in the reaction processes. It is well known that 
the solution of the above equation automatically sat- 
isfies the unitary conditions for the scattering ampli- 
tudes. 

The MB M'B' transition potential is defined 

as 



VmB,M'B'{P,P';E) = VMB,M'B'iP,P') + 

Tn* mb{p)^ n'm' b>{p') 



JV* 



E-m%. 



(2) 



where m^. and Fj^. mb{p) represent the bare mass 
of the «th N* state and the bare N* MB decay 
vertex, respectively. The meson exchange potential 
Vmb,m'b' is derived using the unitary transformation 
method'^'- from the phenomenological Fagrangians 
which respects gauge and chiral symmetries. It is 
noted that the derived potential is independent of 
the total scattering energy E. The second term de- 
scribes a MB M'B' transition through the prop- 
agation of the bare N* state, MB ^ TV* ^ M'B'. 
Also, defining E^ik) — [ni^+k^Y^^ with rria being the 
mass of particle a, the meson-baryon propagators in 
Eq. (dD are GMB{k,E) = 1/[E - EM{k) - EBik) +i£] 
for the stable ttN and r]N channels, and GMB{k,E) — 



1/[E- Euik) - EB{k) - Y.MBik,E)] for the unsta- 
ble ttA, pN, and aN channels. The self energy 
Y,MB{k,E) is calculated from a vertex function defin- 
ing the decay of the considered unstable particle in 
the presence of a spectator tt or iV with momentum 
k. For example, we have for the ttA state. 



mA 



Ea{p) Jo 



q^dq 



M^N{q) 



I/a- 



mAq)+p'V^ 

{q)f 



E-EM-[Ml^{q)+p^]y^ 



, (3) 



where M„ff{q) — E„{q)+Ej^{q) and fA-.-^N{q) defines 
the decay of the A ttN in the rest frame of A. The 
self-energies for pN and aN channels are similar. 

We can split the full partial wave amplitude 
Tmb,m'b'{PtP',E) into two pieces without introduc- 
ing any approximation: 

TmB,M'B'{P,P]E) = tMB.M'B'{P,p']E) + 

tf;B,M'B'{P.P'\E). (4) 

Here the first term is a solution of Eq. ([T|) but re- 
placing the full transition potential Vmb,m' b' with 
the meson-exchange potential vmb.m'b'- Therefore 
tMB,M'B'{p,p'',E) expresses the pure meson-exchange 
processes and thus is called the meson-exchange am- 
plitude in this paper. The second term describes re- 
action processes associated with the bare N* states, 

[D{E)],,,Tr,,,M'B'{p';E), (5) 

where f iv*.jv/'s'(fc;£') is the dressed N* ^ M'B' ver- 
tex function which is defined as 

fiv;,Af's'(p,£') = Tnj,m'b'{p) + 

/ dqq^Tt^.^M"B"{q)GM"B"{q;E)x 

tM"B",Ai'B'{q,p;E). (6) 

The inverse of the propagator of dressed N* states in 
Eq. JS]) is 

[D-\E)],., = (i?-m^.)<5,,-[M(i?)],,„ (7) 
where the N* self-energy is defined by 



[MiE)], 



E 

MB 



q^dqTj, 



Aq.E)y. 



GMB{q,E) {q,E). (8) 



As for the inclusion of the bare N* states, we im- 
pose a condition that the number of the bare N* 
states included should be as minimum as possible. 
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This means that if we can explain the data even with- 
out some of the bare N* states, then we wiU take 
those states out of our framework. On the other 
hand, if we cannot explain some data set with the 
current framework, then we will try to introduce an 
additional bare N* state and it may correspond to 
the new N* state. At present we have introduced 16 
number of the bare N* states. 

As clearly seen in Eqs. ©-([S]), the bare N* 
states couple to the meson-baryon continuum through 
the reaction processes and become resonance states. 
On the other hand, the meson-exchange amplitude, 
which describes the pure meson-exchange processes, 
can also generate the resonance poles dynamically. 
Our model allows both possibilities. 

The MB M'B' amplitude defined in Eqs. U])- 
([Hl constitutes a basic ingredient to construct all sin- 
gle and double meson production reactions with the 
initial ttN and ^'^'^N channels. The details of con- 
structing the TriV — > nnN and ^^"^N MB,7T7tN 
amplitudes can be found in Refs. [3,[H|- 

3 Current status of the dynamical 
coupled-channels analysis at EBAC 

3.1 ttA^ ^ vrA^ scattering 

As steps toward extracting information on the N* 
states, it is necessary to first determine the hadronic 



parameters of the EBAC-DCC model. This was ac- 
complished in Ref. 0|. There the SAID nN partial 
wave amplitudes are considered as the "experimental 
data" and the hadronic parameters are determined 
by fitting to them up to W = 2 GeV. The numerical 
fit was performed systematically using the MINUIT 
library. We refer to this analysis and the obtained 
parameter set as JLMS in the following. 

In Fig. [1] our results of the angular distributions 
and the polarizations of ir^p ■n^p and vr^p 
7r°n calculated from the EBAC-DCC model with the 
JLMS parameters (red solid curves) are compared 
with the SAID results (blue dashed curves). Our re- 
sults agree with the SAID results as well as the data. 
As a next step we plan to improve our hadronic pa- 
rameters by fitting to the data of these observables. 

In Fig. m we also present the resulting total cross 
sections of ttN — > X (red solid curves) and TriV ttN 
(blue dashed curves) which agree with the data well. 
The deviation of the ttN X total cross sections 
above = 1.8 GeV can be understood because in 
those energy region the multi-pion production reac- 
tions, which are not considered in the current EBAC- 
DCC model, also start to visibly contribute to the 
total cross sections. 



jt p— ► n p 
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Fig. 1. DifTerential cross sections (left) and polarizations (right) of ■K~p^-K~p and 7r~j3 — > Tr'^n. The results 
from EBAC-DCC analysis [3| (red solid curves) are compared with those from GWU group Q (blue dashed 
curves) . 
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Fig. 2. The predicted nN total cross sections from the EBAC-DCC model. (Left panel) The n^p-^X (solid 
red curve) and ■K'^p^n'^p (dashed blue curve) reactions. (Right panel) The n'p^X (solid red curve) and 
7r~p— >7r~p + 7r^n (dashed blue curve) reactions. The corresponding data to to the solid and dashed curves 
are presented as open square and open circle, respectively. The data are from Refs. [1, 



3.2 ttA^ ^ vrvrA^ reaction 

With the JLMS parameters obtained from the 
analysis of the ttN scattering, we have calculated 
cross sections of the nN ttttN reactions in Ref . Q . 
Within our model there appears no additional pa- 
rameters in constructing the t:N — > TTTriV amplitude. 
Thus the results are purely our predictions. 

Figure. [3] is the resulting total cross sections (red 
solid curves). Without any modifications of the pa- 
rameters, our results already capture the essential fea- 
tures of the total cross sections up toW — 2 GeV. One 



possible reason for the deviation between our results 
and the data will be because the ttN elastic scatter- 
ing does not completely fix the parameters associated 
with the inelastic ttA, pN , and it TV channels. This re- 
sult indicates that we need a combined, simultaneous 
analysis of the ttN elastic and ttN — > ttttA^ reactions 
to construct more complete hadronic amplitudes up 
toW — 2 GeV. In the same figure, we also present the 
results in which the coupled-channels effect is turned 
off. We observe that the coupled-channels effect is 
significant in all W region. 



+ + + 
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Fig. 3. The predicted ttN nnN total cross sections. Red solid curves are the full results, and blue dashed 
curves are the results in which the coupled-channels effect is turned off. See Ref. [sj for the data references. 
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In Fig. m we present the invariant mass distri- 
butions of TT^p — > TT+TT'n. The shape of our results 
agrees with the data very weU except for the tt+tt" 
distribution at low W. It is expected that this dis- 
agreement with the data provides useful information 
to refine the parameters associated with the coupling 
of the low lying N* states to the pN and aN channels. 



W= 1.44 (GcV) 



W= 1.79 (GcV) 




0.2 0.4 0.6 0.8 
M(jlV) (GcV) 



0.2 0.4 0.6 0.8 
M(7lV) (GcV) 



0.2 0.3 0.4 0.5 0.6 
M(llV) (GcV) 

Fig. 4. The predicted invariant mass distri- 
butions of the TT~p -K^-K'n reaction at 
W = 1.44, 1.66, and 1.79 GeV. The red solid 
curves are the full results, the dotted curves 
are the phase-space normalized to the data. 
The data are from Ref. [l3|, whose magni- 
tude is determined by normalizing them to the 
7r~p — > 7r+7r~n total cross sections listed in 
Ref. 

3.3 jN—^ttN and eN e'nN reactions 

The main purpose for analyzing single pion photo- 
and electro-production reactions is to extract the A^- 
N* electromagnetic transition form factors. Their 
precise determination is crucial for understanding the 
N* structure because their dependence is expected 
to strongly reflect nature of the N* structure. 



Within the EBAC-DCC model the N-N* electro- 
magnetic transition form factors are obtained from 
the dressed 'j*N ^N* vertex function defined by 



E 

M"B'' 



dkk'^v^(,)ff^M„g„{q,k,Q^) x 

Gm"B" {k', E)^ M" B" ,N* 

{k;E), (9) 

where r^t.j^jy. (g, Q^) is a bare 7**'iV — > N* ver- 
tex function and v^(,)j^ j^j,,g,,(k,q,Q^) is a meson- 
exchange j'-*'' N ^ M" B" transition potential. 

Our first analysis of the jp — > ttN and p{e, e'Ti)N 
reactions has been performed up to < 1.6 GeV in 
Refs. 12 1 and [l3|, respectively. In the analysis we 
fixed the hadronic parameters with the JLMS values 
and varied only the parameters associated with elec- 

.{qM^) and 



tromagnetic interactions, i.e. 

V^i.)N,MB{k,q,Q'^)- 

Figures E] and E] present our results of the total 
cross sections of 7p ttN and the five-fold differen- 
tial cross sections oi p{e,e'7T)N at — OA (GeV/c)^, 
respectively (red solid curve). Our results agree with 
the data very well in the considered energy region up 
to Vl^ = 1.6 GeV. We also present that the results in 
which the coupled-channels effect on the electromag- 
netic interactions is turned off (blue dashed curves). 
As for the photoproduction reactions, we find that 
the coupled-channels effect has about 30-40% of con- 
tributions to the cross sections in all W region up 
to 1.6 GeV. On the other hand, the coupled-channels 
effect on the electroproduction reactions is still large 
around W — 1.2 GeV but becomes small rapidly at 
high W for increasing Q^. 

The resulting Gl^, G^, and GJ, form factors of 
7*iV A(1232) and the helicity amplitudes for the 
higher N* states evaluated at their Breit-Wigner 
masses are found in Refs. ij] and 13 1. 



400 



300 




1100 1200 1300 1400 1500 1600 " 1100 1200 1300 1400 1500 1600 
W (MeV) W (MeV) 

Fig. 5. The total cross sections of jp yr^p (left) and jp iv'^n (right). The red solid curves are the full 
results of EBAC-DCC model; the blue dashed curves are the results in which the coupled-channels effect in 
the electromagnetic interactions is turned off (see Ref. [l^ for the details). 
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Fig. 6. The five-fold differential cross sections of p(e,e'7r'')p (upper panels) and p{e,e'-K^)n (lower panels) at 
= 0.4 (GeV/c)^. The meaning of each curve is same as in Fig. [S] The results are taken from Ref. 13]. 



Extraction of A^* poles from the con- 
structed scattering amplitudes 



Once the scattering amplitude is constructed from 
analyzing the data of meson production reactions, we 
can extract various information on the N* states such 
as masses, widths, and decay vertex functions. The 
general scattering theory tells us that the N* mass 
and width should be identified with a pole position of 
the scattering amplitude on the complex energy plane 
and the N* decay vertex functions with the residues 
at the corresponding pole'^'— 

However, the scattering amplitude, which is ob- 
tained as a solution of the coupled-channels equa- 
tion ([1]), is originally not defined on the complex en- 
ergy region where the resonance poles exist. To ex- 
plore the N* pole positions, we need to make an ana- 
lytic continuation of the scattering amplitudes. How 
to perform such analytic continuation both mathe- 
matically and numerically is described in detail in 
Ref. 17[ . Therefore here we just present the results of 
the N* pole positions extracted from our constructed 
scattering amplitude. 

In Table [I] the extracted N* pole positions from 
our current EBAC-DCC model arc compared with 
the PDG values. It is found that most of our results 
and the PDG values are close to each other. It is 
noted, however, that we obtain these pole positions 
from just analyzing the TriV scattering and did not 
try to make our results closed to the PDG values at 
all. 

Within our current model, all the extracted N* 



states are found to be originated from the bare N* 
states [the second term of Eq. (0)] and there finds no 
meson-baryon molecular type of resonances, which is 
generated dynamically from the pure meson-exchange 
processes. 

Table 1. The resonance pole positions m_B 
[listed as (Re mij,— Im mij)] extracted from 
the EBAC-DCC model with the JLMS pa- 
rameter set are compared with the values of 
3- and 4-stars nucleon resonances listed in 
the PDG [i]. "— " for P33(1600), P13 and 
P31 indicates that no resonance pole has been 
found in the considered complex energy re- 
gion, Re{E) < 2000 MeV and -Im(£') < 250 
MeV. 
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We also find in Table [T] that tliere exist two Pn 
poles, (1357,76) and (1364,105), in the energy region 
close to the Roper resonance. Here it is noted that 
recent analysis by the GWU and Jiilich groups have 
also reported the two Pn poles in the same energy 
region-^' . This is quite remarkable because, al- 
though the three approaches are based on completely 
different pictures for the Roper resonance, all of the 
three show the two poles of the Roper resonance in the 
same complex energy region after fixing their model 



parameters by fitting to the ttN scattering data. To 
conclude whether both of the two poles are the con- 
stituents of the Roper resonance, however, we need 
to make clear how these two poles affect the physics 
on the real energy axis. This will be accomplished by 
investigating the magnitude of the residue at the pole 
positions. 

In Fig. [71 we present the 3D plots of the N* pole 
positions on the complex energy plane. 



'11 




D 
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1 400 J 





as 





2000 



1800 ' 1800 ' 1800 

2000 2000 

ReE(MeV) ReE(MeV) He E (MeV) 

Fig. 7. The 3D plots of the N* pole positions on the complex E plane. The first two P\\ resonances are very 
close and cannot be shown explicitly. 

5 Plans for future developments 

As a future work, the EBAC program has four 
main components described in the following subsec- 
tions. 



5.1 Extracting N* 
tex functions 



■MB and N* 



'1 



ver- 



In Sec. U we have presented the N* pole posi- 
tions extracted from the current EBAC-DCC model. 
As a next step we will move to the extraction of the 
N* MB and N* — > j'-*''N vertex functions. As 
mentioned in Sec.[4l a proper, physically meaningful 
N* vertex function is defined as a residue at the cor- 
responding N* pole, not at the Breit-Wigner mass. 



Therefore we need to perform an analytic continua- 
tion of the scattering amplitude as well as in the case 
of pole extraction. 

5.2 Analysis of photoproduction reactions 
and search for new N* states 

We plan to analyze photoproduction reactions of 
various final states such as KY, rjN, and ttttA^ in the 
same way as made in the single pion photoproduc- 
tion reactions. As for these final states, the amount 
of currently available data of the initial 7A'^ are much 
larger than those of the initial ttA^. Therefore we 
will use these photoproduction data to precisely de- 
termine the parameters associated with TV* KY, 
N* rjN , and N* TTA,pN,aN decay vertex func- 
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tions. Also, the analysis of these reactions is interest- 
ing because recent experiments have suggested new 
N* states which strongly couple to these reaction 
channels but not to the dominant ttN channel'^' 
We could find a possible existence of such new N* 
states from the analysis. 

5.3 Developing connection with hadron 
structure calculations 

All N* parameters in the EBAC-DCC model have 
hitherto been determined by fitting to the data. 
While this has greatly improved our knowledge of 
the spectrum, it has not improved our understand- 
ing within QCD. For this, it is necessary to develop 
a method to connect our extracted N* information 
to the various hadron structure calculations. This is 
crucial for reaching the final goal of EBAC. Encour- 
agingly, in that connection progress has been made 
in constituent quark models'^^' and Lattice QCD 
simulations'^'^ . 

5.4 Upgrading EBAC-DCC model 

In parallel with the above investigations, we keep 
upgrading EBAC-DCC model. At present the param- 
eters associated with the hadronic interactions are 
determined by just analyzing the ttA^ ttN scat- 
tering, and the electromagnetic parameters are de- 
termined from the analysis of the single pion photo- 
and electro-production reactions in which the fixed 
value of the hadronic parameters are used. However, 
we have seen that we need a combined, simultane- 
ous analysis of the ttN and ttttN channels to obtain 



more reliable model. In the next step we will refine 
our model by performing the combined analysis of the 
ttN ^ ttN^ttttN and jN ^ ttNjTtttN reactions. 




Fig. 8. One example of diagrams including the 
three-body nnN unitary cut. The intermedi- 
ate -kttN state at the red dotted line can be 
on the mass-shell. 

Also, we will incorporate all processes with the 
TTTriV 3-body unitary cut. Some of the processes as 
shown in Fig. [S] have not been included in our cur- 
rent model. Such processes could be important for 
the meson production reactions up to 2 GeV. As far 
as we know, such complete treatment of the three- 
body unitary cut has never been done before in the 
N* study, and thus this improvement will make our 
model unique from other dynamical models. 
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